arXiv:1507.04889vl [astro-ph.GA] 17Jul2015 


Draft version July 20, 2015 

Preprint typeset using D-T^X style emulateapj v. 5/2/11 


A HYPER SUPRIME-CAM VIEW OF 
THE INTERACTING GALAXIES OF THE M81 GROUP* 

Sakurako Okamoto 1 , Nobuo Arimoto 2,3 , Annette M.N. Ferguson 4 , Edouard J. Bernard 4 , Mike J. Irwin 5 , 

Yoshihiko Yamada 6 , and Yousuke Utsumi' 

1 Shanghai Astronomical Observatory, Nandan road, Shanghai, 200030, China 
2 Subaru Telescope, National Astronomical Observatory of Japan, 650 North A’ohoku Place Hilo, HI 96720, U.S.A. 

3 The Graduate University for Advanced Studies, Osawa 2-21-1, Mitaka, Tokyo, 181-8588, Japan 
institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ U.K. 
institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, U.K. 

®National Astronomical Observatory of Japan, Osawa 2-21-1, Mitaka, Tokyo, 181-8588, Japan 
7 Hiroshima Astrophysical Science Center, Hiroshima University, Kagamiyama 1-3-1, Higashi-Hiroshima, Hiroshima 739-8526, Japan 

Draft version July 20, 2015 

ABSTRACT 

We present the first results of a wide-field mapping survey of the M81 group conducted with Hyper 
Suprime-Cam on the Subaru Telescope. Our deep photometry reaches ~ 2 magnitudes below the tip 
of the red giant branch (RGB) and reveals the spatial distribution of both old and young stars over an 
area of ~ 100 x 115 kpc at the distance of M81. The young stars (~ 30 — 160 Myr old) closely follow 
the neutral hydrogen distribution and can be found in a stellar stream between M81 and NGC 3077 
and in numerous outlying stellar associations, including the known concentrations of Arp’s Loop, 
HolmberglX, an arc in the halo of M82, BK3N, and the Garland. Many of these groupings do not 
have counterparts in the RGB maps, suggesting they may be genuinely young systems. Our survey 
also reveals for the first time the very extended (> 2 x R 25 ) halos of RGB stars around M81, M82 and 
NGC 3077, as well as faint tidal streams that link these systems. The halos of M82 and NGC 3077 
exhibit highly disturbed morphologies, presumably a consequence of the recent gravitational encounter 
and their ongoing disruption. While the halos of M81, NGC 3077 and the inner halo of M82 have the 
similar (g — i )0 colors, the outer halo of M82 is significantly bluer indicating it is more metal-poor. 
Remarkably, our deep panoramic view of the M81 group demonstrates that the complexity long-known 
to be present in HI is equally matched in the low surface brightness stellar component. 

Subject headings: galaxies: groups: individual (M81) — galaxies: halos — galaxies: individual(M81, 

M82, NGC 3077) — galaxies: interactions — galaxies: stellar content — galaxies: 
structure 


1. INTRODUCTION 

Over the last decade, deep studies of nearby galax¬ 
ies have led to the discovery of vas t ste lla r en vel opes 
that are often rich in substr ucture (iMihos et all 120051 
iMartm cz-Dclg ado et all I2010D . These components are 
naturally predicted in models of hierarchical galaxy as¬ 
sembly and their observed properties place important 
constraints on the amount, nature and history of satellite 
accre tion (e.g. iBullock fc .Iohnstonll2005l iPillenich et al.1 
I2014T) . Due to their very low surface brightness, one of 
the most effective ways of mapping the peripheral re¬ 
gions of galaxies is through resolved star studies. For 
example, dedicated surveys of red giant branch (RGB) 
stars around M31 have revealed a stellar halo extending 
to more than ~ 200 kpc that is dominated by tidal de¬ 
bris features (e.g.llbata et alJl200lHFerguson et alJl2002t 
iMcConnachie et all 120091 llbata et al.l 120141) . Similarity, 
the Sloan Digital Sky Survey (SDSS) has been used to ex¬ 
plore main-sequence (MS) turn off stars in the halo of the 
Milky Way, leading to many discoveries of new substruc¬ 
tures, satellites and a refine d characterization of h alo and 
thick disk properties (e.g. IBelokurov et al.l 120071) . 

Using wide-held cameras equipped to 8m class tele¬ 
scopes, it has recently become possible to extend 

t Based on data collected at Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan. 


these studies to systems beyond the Local Group (e.g. 


Mouhcine et al.ll2010l : lBarker et al.f2012t|Crnoievi6 etal 


2013). Located at a distance of 3.6 Mpc dFreedman et al 


19941 ), M81, is a prime target for wide-held mapping 


of its resolved stellar content. Spectacular neutral hy¬ 
drogen images have demonstrated the signihcant tidal 
interactions between M81 and its two brightest neigh¬ 
bors, M82 and NGC 3077, which mo delling suggests have 
taken place in the last 300 Myrs (e.g. Ivan der Hulstlll979l : 

lYun et al.ll994IYunlll999HChvnoweth et al.ll2008l) . Deep 

photometry from the Hubble Space Telescope (HST) 
has been used to argue that the outlying HI concentra¬ 
tions Arp’s Loop (AL), and HolmberglX (HoIX) may 
be tidal d warf galaxies formed a s a result of these in¬ 
teractions ([Makarova et all l2002t Ide Mello et al.l 120081 
ISabbi et al.ll2008l) . In NGC 3077. 90% of the atomic hy¬ 
drogen is located eastward of the center, in the tidal arm 
called “the Garland” where young stars have been ob¬ 
served (Karachentsev et all 1 19851 iSakai fe MadorelfeOOll : 
IWeisz et all 120081 b Several other young star concen- 

trati ons have been as s ociated with peaks in the HI 
gas (iDur rell et, a ll 2004 S i mj 4 _djj 200 ^ [2?To 8 bt 

iMouhcine fc Ibatal l2009l :~ lChiboucas et, al.l 120131b how¬ 


ever, the global properties of this population throughout 
the M81 group are still poorly known. 

The old stellar content arou nd M81 has a l so bee n stud¬ 
ied using large telescopes. I Barker et alJ (120091 ) found 
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the evidence for a faint, extended structural compo¬ 
nent beyond the bright optical disk of M81 from wide- 
held images taken by Subaru/Suprime-Cam. They de- 
tected no colour g radien t in this structure out to 44 kpc; 
iMonachesi et, al.l (1201 - 111 used HST po i ntings to extend 
this result to 50 koc. IChiboucas et al .1 (|2013l ) confirmed 
12 new dwarf satellites as members of the M81 group, 
discovered from a 65 deg 2 survey with CFHT/MegaCam. 

In this paper, we present the first results from a deep 
wide-held imaging survey of the M81 group that we 
are conducting with the new prime-focus imager, Hyper 
Suprime-Cam (HSC), on the Subaru Telescope. We re¬ 
port on the analysis of the inner 4 deg 2 area, correspond¬ 
ing to a region spanning 100 x 115 kpc at the distance of 
the galaxy, which reveals the hrst truly panoramic view 
of the low surface brightness stellar component. The ob¬ 
servations and data reduction are described in Section [2j 
Section [3] and 0] present our analysis and results, which 
are discussed and conclusions drawn in Section [5] 

For this paper, we adopt a distance modulus 
for M 81 and associated sy s tem s of (m — M) 0 = 

27.79 dRadburn-Smith et all 1201 lil . position angles for 
M81, M82, and NGC 3077 of 157°, 67.5°, and 55.0° 
east of north, R 25 radii for M81, M82, NGC 3077 
and HoIX of 13. 8b 5 .6', 2.7 , _jmd_2 ; 5f L _ respective ly 
(Ide Vaucouleurs et all 1199 ll iKarachentsev et ahl 12004 ). 

2 . OBSERVATIONS AND DATA REDUCTION 

We observed the central region of the M81 group in 
the g- and i-bands using 4 pointings of Subaru/HSC 
during the nights of January 21 and 22, 2015 (PI: S. 
Okamoto; Proposal ID: S14B-101) with the seeing ranged 
from 0.6"to 0.9". The HSC consists of 104 CCD de¬ 
tectors and provide s a field-of-view of 1 .76 deg 2 with a 
pixel scale of 0.17" (lMivazaki et al.l 120121) . The observa¬ 
tions were obtained as part of a survey to map the M81 
group with 7 HSC pointings. In this paper, we focus on 
the inner 4 deg 2 of o ur sur vey which overlaps the SDSS 
footprint (lYork et al.ll 2000 ll . 

The raw images were processed using the HSC pipeline 
(version 3.2.2), which is based on a software suite be¬ 
ing developed fo r the Large Synoptic Survey Telescop e 
(LSST) project (jlvezic et, all 120081 : 1 Axelrod et, al.ll2010H . 
For the processed images, the DAOPHOT in IRAF was 
used to obtain the PSF photometry of resolved stars 
dStetsonl 119871) . Astrometric and photometric calibra¬ 
tions were done using the SDSS catalog. Artificial star 
tests were performed on some parts of the reduced images 
using the ADDSTAR in DAOPHOT, and indicate that 
our photometry is at least 50% complete to 26 mag in 
both bands, except for the inner regions of galaxies. We 
separate point sources from extended sou rces i n the same 
mann er as for Suprime-Cam images in Ok amoto et al.1 
( 2012 !). Full details of the observations and data reduc¬ 
tion will be presented in a forthcoming paper. 

3. THE COLOR-MAGNITUDE DIAGRAMS 

Figure Q] shows the resulting CMD of roughly 550,000 
point sources found in the whole 4 deg 2 field. The error 
bars represent the photometric errors at (g — i)o = 0 , 
as estimated by the artificia l star tests. The G alacti c 
extinction is taken from iSchlaflv fe Finkbeinerl (|2011h . 
Since the extinction varies across the observed field, we 



(g-0o 


Fig. 1.— Dereddened CMD of stellar objects located within the 
central 4 deg 2 area. The dashed boxes delineate the selection cri¬ 
teria for different stellar populations and are used to construct 
the maps presented in Figure [3] Theoretical Padova isochrones, 
adjusted at (m — M )o = 27.79, are shown for a 12 Gyr old popu¬ 
lation with [M/H] = —2.2, —1.75, —1.3, —0.75 from the left to the 
right (magenta), and for a [M/H] =-0.75 population with ages of 
10,18,32,50,100,160 Myr from top to bottom (blue). The dotted 
lines represent the completeness levels of 50% and 30%. 


apply a reddening correction to each source indiv i dually 
according to its location, assuming a iFitzpatrickl (119991) 
reddening law with Ry=3.1. The central region within 
r = 15' of M82 in the reddening map shows significantly 
higher extinction E(B — V) ~ 0.16 that includes the in¬ 
ternal reddening of M82. Therefore, we replace it with 
E(B-V) = 0.075. 

To aid in understanding the range of stellar popula¬ 
tions present, theor etic al Padova isochrones are over¬ 
laid {Brcssan et al.1 l2012fl . We find that tracks with 
metallicity [M/H] varying from ~0.0 to below —1.0 for 
young stars, and [M/H] < —1.0 for old stars pro¬ 
vide a good description of the data, in agreem ent with 
other studies of smaller regions ( Makarova et al .1 | 2002 t 
Williams et al.ll2009t iDurrell et al. 201(1 Kndritzki et al l 


20121 iBarker et al.ll2009t IMonachesi et al.ll2013h . We use 

[M/H] = —0.75 as the fiducial value for the young pop¬ 
ulation, and overlay the 10-160 Myr old isochrones as 
blue solid lines in Figure [L] For the old population, we 
plot [M/H] = —2.2 to —0.75 isochrones of 12 Gyr old as 
magenta solid lines. 

Figure [T| is mostly populated by old RGB stars lo¬ 
cated at in j b 24 a nd (g — i) 0 ~ 1 . 2 . As discussed by 
IBarker et al.1 (l2009f ). the over-density at in ~ 26 and 
(g—i)o ~ 0.3 (labelled BO for ‘blue objects’) mostly sam¬ 
ples unresolved background sources. On the blue side of 
the foreground Galactic dwarf sequence (at (g — i ) 0 ~ 
0.4), young MS and core Helium burning (cHeB) stars in 
the M81 group are found. We select stars in different evo¬ 
lutionary phases as shown by the dashed boxes in Figure 
[TJ MS, cHeB, Red Supergiant (RSG), Asymptotic Gi¬ 
ant Branch (AGB) and RGB stars. The boundaries were 
adopted to limit the number of foreground/background 
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Fig. 2. — Dereddened CMDs of stellar sources in the disk and halo (see text) of each galaxy, and in the SE-stream. The overplotted 
isochrones are the same as those in Figure [l] 
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NGC3077 (r<1 Omin) 


SE-slream 


contaminants. The young MS box mainly contains stars 
younger than ~ 50 Myr old while the MS+cHeB box 
is occupied by MS and post MS stars of < 100 Myr old. 
On the red bright side, the polygon contains RSGs about 
25-160 Myr old with some contamination from Galactic 
disk stars. Above the RGB tip at iq = 2 4.25, interme¬ 
diate age (~ 0.5-8 Gyr) AGBs are found dMarigo et al.1 
1200811 . From the RGB tip to about 1.5 magnitude below, 
the blue and red RGBs boxes contain stars older than 
1 Gyr. We also examine the spatial distribution of BO 
sources. We note that the completeness of our photom¬ 
etry decreases towards the red, as shown by the dashed 
lines, so our RGB sample may be biased toward the bluer 
side. 

Figure [2] shows the dereddened CMDs of the disk and 
halo of each galaxy, BoIX and the stream between M81 
and NGC 3077 (hereafter the SE-stream). Stars within 
r = 2 x R 25 (i.e. 27.6' or 29 kpc) of M81, and within 10' 
(10.5 kpc) for M82 and NGC3077, within r = R 25 (2.5' 
or 2.6 kpc) of BoIX are shown (see dashed lines in Figure 
[3D; an axis ratio of 0.75 has been used for M82 to take 
into account the varying flattening of the stellar distribu¬ 
tion in the outer regions. We also select a 15' x 6 ' area for 
the SE-stream between M81 and NGC 3077. Young blue 
MS and RSG can be easily seen in the M81 CMD, as well 
as vast numbers of RGBs. In M82 and NGC 3077, RGBs 
are prominent and some MSs exist, but few if any RSGs 
can be seen. We note that the area within R 25 of each 
galaxy could not be resolved due to crowding, so we miss 
the stars of the M81 and M82 disks, and in NGC 3077 
center. In the SE-stream, the most luminous MSs cor¬ 
respond to ~32 Myr old. In BoIX, several MSs, cBeBs 
and RSGs exist. Although we cannot resolve RGBs at 
the innermost (< 1') of BoIX, the number density of 
RGBs in BoIX are comparable to those of other regions 
at t he same distance from M81. Therefore, as discussed 
bv lSabbi et all (I2008H with deeper BST images, most of 
the old component at BoIX may belong to M81 halo. 

4. TBE SPATIAL DISTRIBUTIONS OF YOUNG 
AND OLD COMPONENTS 


Figure ET] show s the spatial distributions of RI gas 
(lYun et allll994f) . MS, cBeB, RSG, AGB, blue and red 
RGB stars, and BOs defined in Figure [I] without cor¬ 
rection for the completeness and contaminants. In the 
upper panels, the solid ellipses represent the R 25 radii of 
the three main galaxies. 

The distribution of the young populations traced by 
the young MS, MS+cBeB and RSG stars agrees ex¬ 
tremely well with that of the RI distribution, except for 
the stream at the northwest of NGC 3077 where few stars 
are seen. The young MSs are mainly concentrated in the 
spiral arms of M81, at the north-west side (hereafter NW- 
arm), AL, BoIX, BK3N, the Garland, and SE-stream. A 
number of small clumps are also seen as they follow HI 
blobs around these systems. The stellar concen tration in 
the NW-arm was identified bv iDavideel (l2008bll and was 
sugges ted to be part of the M81 arm bv iBarker et al.1 
( 2009 1. We confirm that it is connected to the stel¬ 
lar concentration on the north arm. The higher den¬ 
sity regions in the SE-stream have been reported either 


as clumps or as a dwarf galaxy (d0959 + 68 ) in 
ous studies (IDurrell et al.l 2004)Mouhcine & Ibatal 

previ- 

2009; 

Chiboucas et al. 2013j). As Chiboucas et al. (2013 

) dis- 


cussed, these over-densities are clearly parts of a single 
stream. Near NGC 3077, many MSs are found in the 
Garland and up to about 8 kpc to the south and 10 kpc 
to the e ast where HI gas an d dust emission have been ob¬ 
served (IWalter et alJl20lTh . In M82, a prominent stellar 
feature ca n be seen a t (Aa Ad) ~ (0.1, 0.6 ) ident ified as 
an arc bv lSun et al.1 (|2005l l and iDavidgel (|2008aH . MSs 
are also distributed well beyond the R 25 radius up to the 
projected distance of about 16 kpc from M82. 

Stars in the MS+cHeB box have a very similar dis¬ 
tribution to that of the young MSs. The arc at the 
southeast of M82 appears to have a clumpy shape, but it 
is an artificial appearance due to the overlap of plotted 
points (see Figure 0. Interestingly, our maps show that 
there is another young stellar feature on the opposite 
side of M82 that appears to be aligned with the south¬ 
ern arc and may therefore be related. The selection box 
of MS+cHeB stars includes some contamination, mainly 
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Fig. 3.— The spatial distribution of s tars in each evolut ionary phase selected in Figurc[l] Shown from the top left to the bottom right are 
the HI column density map taken from lYun et al.; (1994 1. the spatial distributions of stars in evolutionary phases and blue contaminants. 
The cross marks represent the centers of known M81 group members. The solid lines are R 25 of galaxies, with axis ratios of M82 and 
NGC3077 of b/a = 0.38 and 0.83, and an inclination angle of i = 58° for M81, respectively. The dashed lines outline the regions used for 
the CMDs in Figure^ 


from background blue objects, as can be seen by the 
low-level uniform distribution of sources through out the 
area. In the top-right panel of Figure [3l the distribu¬ 
tion of the RSGs is almost the same as that of MSs and 
cHeBs. However, the fainter substructures - SE-stream, 
M82-arc, some clumps around HoIX and AL, and BK3N 
- can not be seen in this map due to the shorter lifetime 
and the lower number of RSGs compared to MSs. 

In contrast, the older populations (AGB, RGBs) have 
a much smoother distribution than the younger stars. 
Old stars are mainly embedded in the halos which reach 
far beyond the R 2.5 radii of three galaxies. The sizes of 
these RGB halos are considerable, and they may even 
overlap (see Figure [fj. While the M82 halo seems to be 
more extended in the blue RGB map, the halos of M81 
and NGC 3077 are more prominent in the red RGB map, 
suggesting they have a higher mean metallicity. In the 
contour map of red RGBs in Figured a tidal stream be¬ 
tween M81 and M82 can clearly be seen, and the outer re¬ 
gions of M82 and NGC 3077 exhibit an “S-shaped” mor¬ 
phology. The dwarf galaxies IKN, BK5N, and KDG61 
cannot be seen in the maps of young stars, but appear 
as over-densities of old populations, implying they have 
not formed as a result of the recent interaction. 

In the bottom-right panel of Figure [3] we plot the dis¬ 


tribution of BOs at 25 < io < 26 and 0.2 < (g — i )o < 0.4. 
The uniform distribution of these sources supports their 
identification as contaminants, since the star/galaxy sep¬ 
aration in our photometry degrades at magnitude fainter 
than about 24 mag in both g- and i-bands. 

5. DISCUSSION AND SUMMARY 

We find that the young intra-group population in the 
M81 group traces the filamentary structures of the HI 
gas connecting M81, M82, and NGC 3077, confirming the 
results of several smaller field-of-view studies. The left 
panel of Figure [5] shows the spatial distribution of stars 
in young MS and MS+cHeB boxes of Figure|l] which are 
color coded according to i-band magnitude in a transpar¬ 
ent manner, so that colors of overlapping points repre¬ 
sent the average color. Bright stars are mainly located in 
the inner disk of M81, while most of young stars in AL, 
NW-arm, BK3N, Garland, and other debris features are 
fainter than io ~ 24 and have similar luminosity distri¬ 
butions to the SE-stream, implying ages of 30-160 Myr 
old. This suggests that star formation in these tidal fea¬ 
tures was synchronized, and may have stopped about 
30 Myrs ago. The SE-stream is slightly inclined from 
the southeast to the northwest, whereas the tail of RGB 
stars in NGC 3077 at the same location is extended to- 
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Fig. 4.- The isodensity contour map of red RGB stars, featured 
to faint structures up to 20cr above the background level. The 
kernel density is estimated with the bandwidth of 1.2'. The marks 
and lines are the same as in Figure [3j 

ward the northeast-southwest. That fact that these spa¬ 
tial distributions do not exactly match may indicate that 
the SE-stream comes from material torn from M81 while 
the RGB tail is material stripped from NGC 3077. The 
color of HoIX is slightly greener than other clumps and 
streams, meaning it includes some bright (younger) MS 
stars as we see in Figure [2] 

The right panel of Figure [5] shows the color distribution 
of RGBs, which can be interpreted as a rough proxy for 
metallicity. The bluest (index=0) and reddest (index=l) 
colors correspond to [M/H] = —2.3 and —0.75, respec¬ 
tively, assuming 12 Gyr old age. The RGBs between 
the solid and dashed lines of M81, M82 and NGC 3077 
have the similar colors; the medians of the color in¬ 
dexes are 0.41, 0.40, 0.46 respectively, corresponding to 
[M/H] ~ —1.4, —1.4 and —1.3. This M81 halo metallic¬ 
ity is slightly lower than the value of [M/H] ~ —1.1 de¬ 
rived from previous Subaru imagery and t he estimation 
[Fe/H = —1.2 from dee p HST photometry (| Barker et al.1 
2009; Monachesi et alJ[2013f) . This might due to the 
missing metal-rich RGB stars in our photometry, since 
the completeness gets worse at redder colors. The halo 
of M82 has a color gradient from the inner greener area 
to the outer bluer region; the greener part is also ex¬ 
tended toward the northeast and along the direction of 


the young stellar arc. Note that we do not correct for 
the internal extinction of each galaxy, so it is not possi¬ 
ble from this study alone to determine if these are the 
bona fide features of the M82 halo. The tidal stream be¬ 
tween M81 and M82 is also predominantly blue in color, 
indicating that this is material being stripped from M82 
onto M81. 

In Figure [4] and [5] the NGC 3077 halo is extended far 
beyond the R 25 , and has a rhombus shape stretched in 
the east-west direction. In the outermost region, an S- 
shape distortion can be disc erned, which is sim ilar to 
what was found around M33 dMcConnachie et al.ll2010l) . 
The component in the northwest appears to reach a max¬ 
imum projected radius from NGC 3077 of ~ 65 kpc, 
but does not appear to trace the HI distribution in 
this region. The S-shaped structure is typical of an 
interacting dwarf galax y with a larger companion (e.g. 
IPeharrubia et all [20 091. Numerical modelling suggests 
the encounters between NGC 30 77, M 81 and M82 took 
place ~ 200-300 Myr ago i Vmii Il999h . which may not 
leave enough time to restore equilibrium in the NGC 3077 
halo. We will return to the topic of the stellar streams 
in the M81 group in a later paper. 

The close encounters between M81, M82 and 

NGC 3077 induced star formation in tidally stripped gas. 
As a consequence, new stellar concentrations were born 
out of these HI rich clumps, many of which lie far from 
the main bodies of the galaxies. Of these concentra¬ 
tions, only AL appears to have a clear counterpart in the 
RGB map. The presence of an older stellar component 
suggests that this object, like the dwarf galaxies IKN, 
BK5N, and KDG 61, may not have a tidal origin. The 
gravitational interactions between the M81 group galax¬ 
ies have also significantly perturbed their older stellar 
components leading to disturbed halo morphologies and 
giant stellar streams which appear to connect all three 
systems. When combined with our forthcoming HSC ob¬ 
servations of the west side of M81, these data will allow 
us to determine the true extent and nature of the intra¬ 
group debris and map the halos of the M81 group galaxies 
to unprecedented distances. 
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Fig. 5.— Left : The spatial distribution of MS and cHeB stars that are color coded according to the luminosity with transparency. 
Right :The spatial distribution of RGB stars. The color of each point represents the (g — z)o color of star with transparency. The marks 
and lines are the same as in Figure [3j 
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